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Abstract 



We study the non-equilibrium condensation process in a holographic superconduc- 
tor. When the temperature T is smaller than a critical temperature Tc-, there are two 
black hole solutions, the Reissner-Nordstrom-AdS black hole and a black hole with 
a scalar hair. In the boundary theory, they can be regarded as the supercooled nor- 
C^ \ mal phase and the superconducting phase, respectively. We consider perturbations 

on supercooled Reissner-Nordstrom-AdS black holes and study their non-linear time 
evolution to know about physical phenomena associated with rapidly-cooled supercon- 
ductors. We find that, for T < T^ the initial perturbations grow exponentially and, 
eventually, spacetimes approach the hairy black holes. We also clarify how the relax- 
ation process from a far-from-equilibrium state proceeds in the boundary theory by 
observing the time dependence of the superconducting order parameter. Finally, we 
study the time evolution of event and apparent horizons and discuss their correspon- 
dence with the entropy of the boundary theory. Our result gives a first step toward 
the holographic understanding of the non-equilibrium process in superconductors. 
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1 Introduction 

Over the past decade the AdS/CFT correspondence has received increasing attention 
and many studies have been made on it. It has been beheved that the AdS/CFT correspon- 
dence plays a central role in the study of the strongly coupled region of quantum field theory 
because it is simply described by classical gravity theory on AdS spacetimes ^-^. Recently, 
the duality between the superconductor and gravity theory has been proposed @-§[ as a new 
application of the AdS/CFT correspondence, in which the simplest gravity theory is given 
by Einstein-Maxwell-charged scalar theory with negative cosmological constant. In this 
gravity theory one of static solutions is the well-known Reissner-Nordstrom-AdS black hole 
solution, in which the scalar field vanishes. It is known that this solution is unstable when 
the temperature T is lower than a critical temperature T^. In consequence of the instability, 
it is expected that the scalar field will condense into a non-vanishing profile and eventually 
the black hole will have the scalar hair breaking the ?7(l)-gauge symmetry spontaneously. 
For low temperature T < Tc, such a static solution has been constructed numerically. It 
was shown that the solution has similar properties with superconductors |^. Thus, the 
instability of the Reissner-Nordstrom-AdS black hole for low temperature was identified 
with the superconducting phase transition. In the same way, the Reissner-Nordstrom-AdS 
and hairy black holes in the gravitational theory were identified with the normal and super- 
conducting phases of a superconductor realized within the dual field theory, respectively. 
Such holographic superconductors are considered as a hopeful approach to understand the 
property of strongly correlated electron systems. 

The non-equilibrium process of strongly correlated systems such as superconductors is 
not fully understood because of difficulties in its theoretical treatment, and has been at- 
tracting much attention |]7|-pJ}1. The AdS/CFT correspondence offers a novel approach to 
this longstanding problem. To understand non-equilibrium process of strongly correlated 
systems, we should simply solve classical dynamics of gravitational systems in the bulk 
thanks to the duality. Many attempts in this direction have been done to get insights for 



strongly coupled quark-gluon plasma, relativistic hydrodynamics and so on [lT2|-pl[]. For 
near equilibrium dynamics of superconductor, there are several approaches from the holog- 
raphy. For example, transport coefficients, such as the electric conductivity, were obtained 



by the linear response theory p,|^,^. (See also [^,^ and references therein.) The static 



and dynamic critical phenomena were also studied in [^.[] These studies played important 
rolls in understanding the holographic superconductor. However, very little progress has 
been made in the regime that the theory is far from equilibrium. In this paper, we give a 
holographic approach to understand the non-equilibrium process of superconductors. 

As we mentioned before, for low temperature T < T^ there are two phases of black 
holes in the gravity side, the Reissner-Nordstrom-AdS and the hairy black holes, and they 
are regarded as supercooled normal phase and superconducting phases, respectively. In 
this paper, we consider small perturbations on the supercooled Reissner-Nordstrom-AdS 
black holes and study non-linear time evolution of them. Because of the instability of the 
Reissner-Nordstrom-AdS black holes, the initial perturbations will grow exponentially at 
the beginning of the time evolution. It is expected that the exponential growth is saturated 



^ In addition, the superfluid hydrodynamics of this system has been discussed in Efl, 27 



due to the non-linear effect and the spacetimes will approach static solutions. It is believed 
that the final states of the time evolution are the hairy black holes obtained in p, but 
there is no proof.0 Studying non-linear time evolutions of the Einstein-Maxwell-charged 
scalar system, we will show that the final states are given by the hairy black holes. We 
are also interested in dynamics of the phase transition on the boundary theory. To reveal 
the non-equilibrium process, we observe the time dependence of the superconducting order 
parameter and study how the normal phase goes to superconducting phase in the boundary 
theory in the middle of the time evolution. We also evaluate the relaxation time scale of 
the boundary theory. Finally, we study the time evolution of event and apparent horizons 
and discuss the correspondence with the entropy of the boundary theory. 

The organization of this paper is as follows. In Section H, we introduce the gravity 
theory of the holographic superconductor and give equations of motion. Assuming the 
plane symmetry, we obtain the (l-l-l)-dimensional partial differential equations (PDEs). 
We also give the outline of our calculations. In Section ^ we solve the equations of motion 
near the AdS boundary in order to obtain boundary conditions at the AdS boundary. In 
Section ^ we explain the numerical method to solve the equations of motion. In Section ^, 
we give the numerical results of non-linear time evolution. We find that the final states 
of the time evolutions coincide with the hairy black holes for general initial conditions. 
In the middle of time evolution, we study the time dependence of superconducting order 
parameters and measure the relaxation time of them. The time evolutions of event and 
apparent horizons are also studied. The final section is devoted to conclusions. 



2 Equations of motion 

We consider the 4-dimensional Einstein-Maxwell-charged scalar theory with negative 
cosmological constant, whose action is given by 



S = d'^Xy/^ 



R+^- -/.'^F'" - I^M^ - ^^^M^l' - ^'l^P 



(2.1) 



where L is the AdS curvature scale and q is the ?7(l)-charge of the complex scalar field z/'- 
The field strength is defined by F = dA as usual. This theory is introduced in 0-§[ as a 
gravity dual of a superconductor. In addition to the diffeomorphism symmetry, this action 
has the local U{1) symmetry, 

A^A + dX, ip^ e^«V , (2.2) 

where A is an arbitrary scalar function. Hereafter, we take the unit of L = 1. 

For simplicity, we assume that the spacetime has the plane symmetry. Using the diffeo- 
morphism and f/(l) gauge symmetries together with the assumed plane symmetry, we can 
take the metric ansatz without loss of generality as 

ds"^ = — - [F{t, z)dt^ + 2dtdz\ + $(t, zfi^dx^ + dy"^) , 

A = a{t,z)dt, (2.3) 

ip = ilj{t,z) . 



■Recently, evidence of this conjecture has been given in 28 1 in the near critical temperature regime. 



We use the ingoing Eddington-Finkelstein coordinates, {t,z). In these coordinates, the 
AdS boundary is located at 2; = 0. These coordinates are convenient for our numerical 
calculations since we can easily extend time slices defined by constant-t surfaces into the 
inside of the event horizon and also set the AdS boundary to be a constant-2; plane. Note 
that diffeomorphism and U{1) gauge symmetries have not been completely fixed, because 
the form of the variables ( |2.3| ) is invariant under the residual symmetry, 1/z — > 1/z + g(t), 
«—)■« + dt\{t) and ip — )■ e'^'^^^^'^ip. These residual gauge symmetries will be fixed by the 
boundary conditions. The complete set of the equations of motion are given by 

($/}$)' - ^ Qz^a'2 + rn'\^\^ -6^=0, (2.4) 

2z^{D^y + iqz^a'^ + 2z^<^~\D<!>)f + 2z^<^~^<!>'Dip + m^ip = , (2.5) 

{z\z"^F)y - z^a''^ + 4$-2(D$)$' - {iIj*'DiIj + ^' D^*) = , (2.6) 

2z'^{Day + z\z-'^Fya' + 4z'^^~\D^)a' - 2zg(^D^* - tlj*DtP) = , (2.7) 



and 



(l>-^Ci = -2$"iL)2$ - (f' - -F J $"^L)$ - |D^p = , (2.8) 

$-2(72 = -2z^<^~\z<^" + 2$') - z^lip'l"^ = , (2.9) 

$-2(^3 = -z^[z'^a" + 2za' + 2z^^~^^'a' + iqiip^*' - ?/^V)] = , (2.10) 



where ' = dz and derivative operator D is defined as 

D^ = dt^-Fdz^/2 , D'^^ = dt{D<^)-Fdz{D<^)/2 , Da = dta - Fdza/2 
Dip = dtip - Fdzi^/2 - iqaip , Dip* = dtip* - Fdzi^*/2 + iqaip* . 



(2.11) 



Here, the operator dt — Fdz/2 represents the derivative along the outgoing null vector. We 
regard Eqs. ( |2.4| - |2.7| ) as evolution equations and Eqs. ( |2.8| - |2.10| ) as constraint equations. 
The constraint equations satisfy 

dzCi = , (2.12) 

DC2 = z^dz{Flz^)C^ + ]^z\dza)C^ , (2.13) 

DC^ = ^z'dz{F/z')Cs , (2.14) 

where DCi = dtCi — FdzCi/2 (z = 2,3). To derive the above equations, we have used 
evolution equations ( |2.4| - |2.7]) . As we can see from Eqs. ( |2.12[ - pTT^ , the evolution equations 
guarantee that the constraint equations to be satisfied if Ci = and C2 = C3 = are 
satisfied on the AdS boundary [z = 0) and on the initial surface (t = 0), respectively. 
Therefore, we solve the evolution equations for the time evolution and use the constraint 
equations only at the AdS boundary and the initial surface to give the boundary conditions 
for the time evolution. 

One of static solutions of the equations of motion is merely the Reissner-Nordstrom-AdS 
black hole solution (with planar horizon), which is given by 

F = 1- 2Mz^ + -Qh"^ , $ = -, a = Qz , ip = . (2.15) 



Parameters M and Q are proportional to mass and charge of the black hole, respectively. 
For this solution, the complex scalar field has a trivial configuration ■?/; = 0. That is, this 
black hole has no "hair" except for the mass and electric charge. The event horizon is 
located dX z = Zj^ determined by F{z^) = 0. In terms of the horizon radius z_|_, we can 
rewrite M as 

1 ^- 1'.2 4 



M=^{l + -^Q'4] ■ (2-16) 



The temperature of the black hole is given by 

1 dF 



T = 

An dz 



=z+ 



12 - 4g2 

±^ . 2.17 



It is known that this solution is unstable when the temperature T is smaller than a critical 
temperature Tc- The numerical value of the Tc is obtained in @-§l. In consequence of 
the instability, it is expected that the scalar field will grow and eventually the black hole 
will have the scalar hair. Thus, the ?7(l)-gauge symmetry is spontaneously broken due 
to condensation of the complex scalar field. Indeed, for low temperature T < Tc at fixed 
charge Q, static solutions of the hairy black hole with ip ^Q exist and have been constructed 
numerically [Q . Similarly to the previous case, the temperature of the hairy black hole is 
given by 



1 dF 

An dz 



(2.18) 



where F(z>) = 0. 



This instability was identified with the superconducting phase transition in the dual 
theory ^^. In this paper, we will investigate the dynamical process of the superconducting 
phase transition in the view of the gravity theory. We show the schematic of our setting in 
Figure |1]. Since we are taking the ingoing Eddington-Finkelstein coordinates ( p.3|) , our time 
slices are given by null surfaces. We consider a slightly-perturbed Reissner-Nordstrom-AdS 
spacetime as initial data, and study non-linear time evolution from it. At the AdS boundary, 
we give appropriate boundary conditions, as we explain in the subsequent section. Inside 
the event horizon, we excise a region before encountering the singularity for the numerical 
calculation. The detailed procedure will be explained in following sections. 



3 Asymptotic expansion 

Now, we determine the asymptotic form of the variables at the AdS boundary z = to 
clarify the boundary conditions for the time evolution. Hereafter, we set the mass of the 
complex scalar field as m'^ = — 2L~^ following 0,^. Then, we can expand variables as 

F{t, z) = l + Fi{t)z + F2{t)z^ + F3{t)z^ + ■ ■ ■ , 

$(t, Z) = 1/z + $o(t) + <^l{t)z + ^2{t)z^ + $3(^)2;' + • • • , 

a{t, r) = ao{t) + ai{t)z + a2{t)z'^ + as{t)z^ -\ , 

2 , „;. U\^3 



^lj{t, r) = Mt)z + ^Ij2{t)z' + ^3(t)z^ + 



Excision of 
the singularity 




event 
horizon 



initial surface 

Reissner-Nordstrom-AdS 
+ perturbation 



Figure 1: Schematic of our setting. We take the ingoing Eddington-Finkelstein coordinates, 
in which time shces are given by null surfaces. We consider a small perturbation on the 
Reissner-Nordstrom-AdS spacetime as initial data on an initial surface and study its non- 
linear time evolution. 



Using residual gauge symmetries Xj z -^ 1/ z + g{t) and a — > a + dt\{t), we put $o(^) = 
and ao(^) = 0. Substituting Eq. ([Ol ) into Eqs. ( p.4H2.10|) and solving the equations order 
by order, we obtain the expansion coefficients as 



1 



Fi=0, F2 = 


--^iT/^t , ^3=2 [^1^1 + ^1^1* 


$1 = ~\ U , 


$2 = --(V^lV^2*+^lV^2) , 

6 


$3 = -^^2p- 

6 


-^^l|'-^(^l4*+A>2) 



V^lV'2 - ^1^2] 



(3.2) 






^2 



--iq{-tp2i^l + ^iV^2) 



where ' = d/dt. We find that all the coefficients Fi{t), $j(t), aj(t) and ipi{t) can be 
expressed by tpi{t) and ip2{t)- This means that if the two functions tpi{t) and ip2{t) are 
given at the boundary we can determine the bulk in principle. It is noticed that these 
are order parameters characterizing normal and superconducting phases in the boundary 
theory [^,0. In our calculation, we do not give ip2{t) at z = 0, and ip2{t) is determined as a 
result of the time evolution. As for the inner boundary, a constant-2 surface, no boundary 
condition is required there since it will be a spacelike surface as long as we take it behind 
the horizon. 

Note that initial values of F^{t) and ai(t) are not determined by the asymptotic ex- 
pansion, because the asymptotic expansion gives just time derivative of F^it) and ai(t). 
In fact, F^lt) and ai{t) are related to with mass and charge of the system. Hence initial 
values of F^(t) and ai(t) represent the initial mass and charge of black holes. We put the 



initial mass and charge as Mq and Qo, respectively and, then, F3 and ai are written as 
F3 = - 2M{t) 

= -2Mo + ^f dt [r,^, + ^^ri + ^iv^2* + rM - \ Ar2 + ^i>2 

ai =Q{t) 

=Qo -iq dt[- 1/J21PI + 1P11/J2 + V'lV^i - V^iV'i] , 



where M{t) and Q{t) denote the total mass and charge, respectively. 



4 A method to solve equations of motion 



(3.3) 



In this section, we explain a method to solve the equations of motion numerically^. From 
Eq. ( ^.2p , we can see that the first few expansion coefficients in Eq. ( \i.l\ ) do not depend 
on ip2{t) but only on ipi{t). In our calculation, we give ipiit) ^^ ^ boundary condition at 
z = 0. Moreover, some variables have divergent parts whose form is fully determined by the 
asymptotically-AdS boundary conditions. Thus, it is reasonable for numerical calculations 
to define new regular variables $, ip and F as 



^p = ziPiit) + z^iP{t, z) , 

F = 1 + zF,{t) + z^F^it) + z^'Fit.z) 



(4.1) 



where $i(t), Fi{t) and F2{t) are given by Eq. (|3.2| ) and they do not depend on ip2{^) but 
only on ifjiit). Note that by these definitions we have 'ip2{t) = ip(t,z)\z=o- We also define 
D$ and D^ as 



D<^ = 1/(2^2) + $i(t)/2 + F2{t)/2 + zD<l>{t, z) 
Dip = -ipi{t)/2 + zDip{t, z) . 



(4.2) 



From Eq. ( ^.11 ), D$ and Dip can be written in terms of $ and ip as 



2D$ = 2<i>i - 2<i> + F + z{2dt^ - d,^ - F2^i) 

- z^{2F2^ + ^iF) - z^{F2d^^ + 2F<I) - z^Fd^^ 
2DiP = 2iPi -2iP + z{2dtip - d,ip - F21P1) 

- z\2F2iP + iPiF) - z\F2d,iP + 2^F) - z^Fd.iP . 



(4.3) 



(4.4) 



In terms of $, il', F, $, ip, D^ and Dip, we can rewrite evolution equations ( |2.4| - p.7| ) as 



{my + Xii^m + X2[^, ^, a] = , 

(D>)' + F[D$,$,V^,a] = , 
(z^F)" + Z[D<^, Di>, $, ^, a] = , 
(Da)' + W[D^, Dip, $, iP,a,F] = , 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



^Similar calculations were performed in [El| in the context of the AdS/QGP. 

6 



and constraint equations ( p. 91) and ( p.lO| ) as 



<I" + Pi[<l>]l>' + P2[<i',^] = 



a 



Qima' + Q, 



. 



(4.9) 
(4.10) 



where Xi, X2, Y, Z, W Pi, P2, Qi and Q2 are functions of each arguments. Since the 
expression of these functions are tedious, we do not write their exphcit expressions. 

The procedure to solve Eqs. ([4.5| - [4.10]) as follows. To begin with, we must construct 
initial data of $(t = 0,z), il){t = 0,z) and a{t = 0,z) on the initial surface t = that 
satisfy the constraint equations ( [4.9| ) and ( |4.1CI| ). 



1. On the initial surface t = 0, we prepare initial data of ip{t = 0,z), initial mass Mq 
and charge Qo- 



2. Integrating Eq. ( [4.9|) and ( [4.10|) from z = to z = zq, we can obtain initial data 
$(t = 0, z) and a(t = 0, 2). Here, z = zo is the boundary of computational domain, 
which we set to be between the horizon and the black hole singularity. For these 
radial integrations, we need boundary conditions $(t = 0,2)1^=0 and a(t = 0,z)\z=o- 
They are given by the asymptotic expansion (|3.1|) as 



l>(t 
a(t 



0,z) 
0,z) 



\z=0 



2=0 



$2(0), dMt = o,z)\^^^ = <^3{o) 

, dza{t = 0,z)\ = ai(0) , 



(4.11) 



where coefficients $2(0), $3(0) and ai(0) are written by '?/'i(0) and 4'2{0) as shown in 
Eq. (|3.2|) . The V'2(0) can be read off from ip{t = 0, z)\z=o. 

Once we have obtained the data of $(t, z), ip{t, z) and a{t, z) on a constant-t surface at 
an arbitrary time t, we can calculate the data of <l>(t + 6t, z), ip{t + 6t, z) and a(t + 6t, z) 
on the next time slice t + (5t as follows. 



i. Integrating Eqs. ( [4.5[ - [4.8|) from z = {^ to z = zq, we obtain the D$, Dip, Da and F on 
the constant-t surface, on which $(t, z), ilj{t, z) and a{t, z) have been known. Boundary 
conditions for the radial integration are given by the asymptotic expansion (p.l|) as 



-^2{t) + ^F.it) - \{Mt)m) + rimiit)) 






1 . N 



(4.12) 



z^F{t,z) 



\z=0 



0, dz[z^F{t,z)] 



2=0 







where $2(^)5 -^3(^) and ai{t) are written by 4>i{t) and V^2(^) as in Eq. (p. 21) and (|3.3|). 



ii. We calculate dt^(t,z), dta(t,z) and dtipit^z) from Eq. ( |2.11| ). We use the upwind 
differencing scheme for advection terms in Eq. (|2.11|). 



^ We do not use Eq. (2^), one of the constraint equations, hereafter. This equation has been aheady 



used to derive the asymptotic form of the variables (3.2). In our calculations, Eq. (2.8) will be guaranteed 



to be satisfied by the boundary conditions (3.2). 



iii. We obtain variables $(t + 5t, z), tp{t + 6t, z) and a(t + 5t, z) on the next time slice t + 5t 
by using time derivative of each variables dt^{t, z), dta{t, z) and dtip(t, z). 

By repeating the above procedure i ~ iii, we can numerically calculate time evolution of 

the system. 

5 Non-equilibrium Condensation Process 

In this section, we summarize the numerical results obtained by the evolution scheme 
of the previous section. In Section |5.1| , we describe the free parameter and the initial 



conditions. We summarize the general properties of the bulk field dynamics in Section ^ 
After that, we study the dynamics of the order parameter of the boundary theory, which is 
the boundary value of the bulk scalar field, in Section ^.3| . We also investigate its growth 



and decay rates around the initial and final states in Section |5]J. In Section |5.5| , we focus 
on the dynamics of the event and the apparent horizons. If a black hole is stationary, then 
the entropy of the boundary theory, which is another important physical quantity, will be 
straightforwardly identified with the horizon area in the bulk. However, it is not so obvious 
in dynamical cases. In order to argue their relevance to the boundary theory entropy, we 
clarify to what extent the bulk spacetime is dynamical and give some consideration on how 



to identify the boundary and the bulk spacetime in Section 5.6. 



5.1 Initial data and parameters 

As explained in Section ^ we should specify the ^(t = 0,z), M and Q on the initial 
surface. We treat M and Q as fixed parameters in this section since they are conserved 
quantities in our setting as we will see later. Without loss of generality, we can fix one of 
the parameters using the scaling symmetry^, 

(t,z,x,y) ^ {kt,kz,kx,ky) , (5.1) 

F ^ F , $ ^ $/A; , a^a/k , tp ^ ^J , (5.2) 

M -^ M/k^ , Q -^ Q/k'^ , T^T/k . (5.3) 

Using this scaling symmetry, in our numerical calculation, we put 

M = i (l + \q'^ . (5.4) 

This condition implies that the horizon radius of the initial black hole is set to unity (see 
Eq. ( p.l6p )p|. As for the initial data of ^, we consider the Gaussian perturbation on the 



Reissner-Nordstrom-AdS spacetime as 

A r r^ _ r ^2-| 

(5.5) 



A 

^{t = 0,2;)= exp 

y ITX 



y )2 



252 



■'This is just a residual coordinate transformation which preserves the form of metric and gauge field in 



Eq. (2^) as well as the boundary conditions F{t,z) ~ 1 and $(t, z) — l/z for z ^- 0. 



6n 



To be accurate, the z^ defined by Eq. (ET6) slightly differs from the apparent horizon position for our 



initial data, because we will give small perturbations on the Reissner-Nordstrom-AdS solutions. 



with A = 0.01, 6 = 0.05 and Zm = 0.3. We tried several other initial conditions and found 
that they yield qualitatively the same results. Thus, we will show below the results only 
for the initial data given by Eq. 

5.2 Dynamics of bulk fields 

First of all, we show the dynamics of the bulk scalar field. As we mentioned before, 
our numerical calculation is performed under the boundary condition ipi (t) = at the AdS 
boundary 2 = 0. This means that we regard ip2{t) as the order parameter of condensation. 
Furthermore, it turns out from Eq. ( p.3| ) that F3 and cti become time- independent when 
ipi{t) = 0, namely the total mass M and the total charge Q are conserved in our calculation. 

In Figure ^ we depict the dynamics of the amplitude of the complex scalar, \ip{t, z)\, for 
q = 1.0 and T/Tc = 0.5 at the initial state. The critical temperature Tc is evaluated for a 
fixed charge Qf\. We can find that much of the wave packet of the initial perturbation ( p.5[ ) 
is instantaneously absorbed in the horizon within tTc < 0.06. Because of the remnant of 
the wave packet, which is the unstable mode contained in the initial perturbation i f).5\j , the 
scalar density grows exponentially for tT^ < 6. The exponential growth is saturated by the 
nonlinear effect at tTc ~ 6. In tTc ^ 6, the scalar density approaches a static solution. As 
in Figure 0, we can find that the static solution coincides with the hairy black hole solution 
obtained in 0. Thus, our result gives numerical proof of the conjecture that the final state 
of the instability of the Reissner-Nordstrom-AdS black hole is the hairy black hole. Our 
result also implies that, for the plane-symmetric perturbations, the hairy black holes are 
stable. It is worth noting that this phase transition from the initial Reissner-Nordstrom- 
AdS black hole to the final hairy black hole is a dynamical process under the fixed mass 
and charge. It implies that the temperature of the initial state and that of the final state 
are different in general. Indeed, the temperature of the final hairy black hole increases 
compared to the initial temperature due to the phase transition. 

5.3 Dynamics of the order parameter 

In following subsections, we show some results from the numerical solutions that are 
relevant to the dual theory. In this subsection, we describe the non-linear dynamics of the 
order parameter of the boundary theory. 

From the asymptotic form of the numerical solution ilj{t, z), we can read oSip2{t) defined 
by Eq. ( p.l| ). The coefficient ip2{t) is regarded as the superconducting order parameter in 
dual theory. Following 0,0, we define the order parameter on the boundary theory as 

{02{t)) = V2ij2{t) . (5.6) 

In Figure ^, we depict the time dependence of the order parameter {q\ {02{t))\Y^'^/Tc, which 
is invariant under the scaling (|5.3|) , for q = 1.5 and T/Tc = 0.2, 0.4, 0.6, 0.8, 1.1, 1.2 and 



^In this section, we investigate the time evohition of the scale invariant variables under the scaling 



symmetry (5.3). For q = 1.0, 1.5 and 2.0, the critical temperature Tc is respectively given by Tc/y/Q 
0.03589, 0.08421 and 0.1234. 





0.08 



0.06 
0.0^ 
0.02 <^ 



(a) < tTc < 14 



(b) 0<tTc< 0.08 



Figure 2: The dynamics of the scalar field for g = 1.0 and initial temperature T/T^ = 0.5. 
In Figure (a), we depict the dynamics of the amplitude of the complex scalar field, \ip(t, z)\, 
on (t, 2;)-plane for < tTc < 14. Because of the instability of the Reissner-Nordstrom- 
AdS black hole, the scalar density grows exponentially for tT^ < 6. We find that, for 
tTc ^ 6, the scalar density approaches some static function. In Figure (b), we depict 
|'(/'(t,2)| for < tTc < 0.08 in order to focus on the behavior of the wave packet of the 
initial perturbation. We can see that the wave packet is refiected by the AdS boundary 
at t ~ 0.04 and much of the wave packet is absorbed in the black hole horizon within 
tTc < 0.06. 



o 




Figure 3: The function r^lipl"^ / Q"^ is depicted on fixed time slices, where r is circumference 
radius defined by r = $. From bottom to top, the curves correspond to tTc = 5.15, 6.44, 
7.73 and 9.02. The top solid curve correspond to that of the hairy black hole. We can see 
that the solution approaches the hairy black hole. 
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1.4. The sharp signal at the small t is caused by the initial perturbation (|5.5|) . For T > Tc, 
the initial perturbation dissipates and the order parameter converges to zero. On the other 
hand, for T < Tc, the order parameter grows exponentially and approaches a non-trivial 
value. We find that, for T < T^ the more rapidly the order parameter converges to its final 
value for lower temperature. 





Figure 4: The dynamics of the order parameter is depicted for g = 1.5. In Figure (a), the 
curves from top to bottom correspond to T/Tc = 0.2, 0.4, 0.6 and 0.8. In Figure (b), the 
curves from top to bottom correspond to T/T^ = 1.1, 1.2 and 1.4. Note that the vertical 
axis in Figure (b) is logarithmic scale, while that is linear scale in Figure (a). 



5.4 Growth and decay rates 

In this subsection, we estimate the growth and decay rates of the order parameter. At 
the beginning of the time evolution, namely, around the normal phase {02{t)) = 0, we fit 
the time dependence of the order parameter as 



|(C»2(t))|=Cexp(-t/trelax) 



(5.7) 



where C and treiax are constants and treiax represents the relaxation time scale. In Fig- 
5(a)| , we depict the 1/treiax around the normal phase against the temperature for q 



ure 



1.0, 

1.5 and 2.0. For T < T^, 1/treiax has negative value since the initial Reissner-Nordstrom- 
AdS black holes are unstable. At the onset of the instability T = Tc, the 1/trciax becomes 
zero. It seems that, for T/Tc — )■ 0, l/treiax dose not diverge, but converges to finite value. 

For T < Tc, the order parameter condenses into the non-trivial value. Around the 
condensed phase, we fit the time dependence of the order parameter as 



|(02(t))|=Ciexp(-t/t,.elax)+C2 



(5. 



where Ci and C2 are constants. In Figure |5(b)| , we depict the 1/trciax around the condensed 
phase against the temperature. Note that we used in this figure the temperature of the 
final state of the time evolution, that is, the hairy black hole. We find that, at T — )■ Tc, 
1/^rciax approaches zero. On the other hand, for T/Tc -^ 0, l/trdax has larger value as the 
temperature becomes lower. 
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In the gravity side, the 1/trciax is nothing but the imaginary part of the quasinormal fre- 
quency of the fundamental mode. For the normal phase, which is the Reissner-Nordstrom- 
AdS black hole, quasinormal modes are studied in [^,^,0. For the condensed phase. 



which is the hairy black hole, quasinormal modes are studied in the decoupling limit g — ?■ oo 
22| . The result in Figure |5(b)| gives the quasinormal frequency of the hairy black holes 



m 



with the back-reaction on the gravitational field. 





0.5 0.6 0.7 0, 



0.9 1 
T/T, 

(a) Normal phase 



1.1 1.2 1.3 1.4 



0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 

T/T, 

(b) Superconducting phase 



Figure 5: The relaxation times treiax of the order parameter against the temperature T. 
Figure (a) and (b) shows the relaxation time around the normal and the condensed phases, 
respectively. The solid, dashed and doted curves correspond to g = 1.0, 1.5 and 2.0. In 
Figure (b), we used the temperature of the final state of the time evolution, that is, the 
hairy black hole. 



5.5 Evolution of horizons 

Now, we investigate the time evolution of apparent and event horizons. The apparent 
horizon z = 2ah(^) can be determined from 



D^t, ZAu{t)) = . 



(5.9) 



We determine the event horizon z = ZEuit) as follows. For sufficiently late time, spacetimes 
settle static solutions. Thus, at late time, the event horizon can be easily determined by 
F{t,ZE}i{t)) = 0. To determine the event horizon for any t, we solve the null geodesic 
equations. 



ZEu{t) 



'-F{t,ZEu{t)) 



(5.10) 



backward along the tangent to the event horizon at late time. Then, we obtain null geodesic 
generators of the event horizon and find the location of the event horizon z = ZEuit). Once 
we know the 2eh(^) and zp^nit), we can calculate the area of event and apparent horizons 

as 



Area(event/apparent horizon) = $(t, 2;eh/ah(^))^ 



(5.11) 
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In Figure P, we depict the time evolution of the area of the horizons. We find that the 
area of the horizons monotonically increase by the time evolution and the event horizon 
has larger area than that of the apparent horizon. We also see that, after the final state has 
settled to equilibrium, the apparent horizon and the event horizon coincide. In addition, 
the two horizons seem to coincide even at the initial state. Thus, it is quite likely that we 
can regard the area of the horizon as the entropy at the initial state. This issue will be 
examined further in the following subsection. 



50.5 



C3 
03 




Figure 6: The time evolution of the area of event and apparent horizons are depicted for 
g = 1.5 and T /Tc = 0.4 at the initial state. Solid and dashed curves correspond to the area 
of event and apparent horizons respectively. 



5.6 Bulk dynamics and bulk-boundary identification 

Having shown properties of the horizons in the bulk, we will discuss how to identify the 
bulk with the boundary based on observations of the bulk field dynamics. When the system 
is stationary, we know that the area of the black hole horizon in the bulk corresponds to the 
entropy of the dual boundary theory. However, there is subtlety in this correspondence for 
dynamical systems (see, for example, [pl|-|33| and so on). When the system is dynamical, 
we do not know the mapping between horizons and AdS boundary a priori. In other 
words, there exists ambiguity of the time slices with which we identify the horizon and the 
boundary. The null time slices we used in this paper, which are illustrated in Figure |l], 
is nothing but one of the uncountable set of time slices. Thus, we cannot immediately 
assert that Figure ^ is equivalent to the dynamics of the entropy in the dual theory. In this 
subsection, we would like to discuss this issue of the identification 

Figure shows time dependence of a metric function <l>. In this figure, we show the value 
of ^'^ /{zT^Y in gray-scale. The dark region, in which (^"^ /{zT^Y is large, represents non- 
stationary region. Rigorously speaking, this figure indicates to what extent a coordinate 
vector ^ = d/dt satisfies the Killing equation: C^g^y = V^C,^ + V^C.^! = 0. The quantity 
IC^Qfj^i^l^ = 8$^/$^ is, however, not always a good indicator of non-equilibrium property 
since ^ coincides with the timelike Killing vector on the AdS boundary and |£^(7^,yp tends 
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to zero there. In order to examine non-staticity even at the boundary, we focus on the 
value of (^/z instead of $ (see Eq. ( p.2| )). As we mentioned before, the initial state and 
the final state are close to equilibrium. It turns out that this fact appears also in Figure ^ 
Another point we should note is that the shaded region extends almost vertically in this 
figure. Since we are using null time-slicing, the above fact means that the non-stationary 
region extends along the ingoing null direction from the boundary to the horizon in the 
bulk spacetime. Consequently, it may be reasonable to suppose that the mapping between 
quantities of the boundary and the horizon is given by the nuU time slices in the current 
case. It will be interesting to examine whether we can apply this idea to other dynamical 
systems. 




Figure 7: Time dependence of the bulk spacetime for the same parameters used in Figure ^ 
The value of ^"^ /{zT^Y is plotted in gray-scale, in which the dark part represents the 
non-stationary region. In order to examine non-staticity of the entire bulk including the 
near-boundary region, we have focused on ^/ z, which takes finite value even on the AdS 
boundary, rather than on $. The dashed curve represents the apparent horizon. 



6 Conclusions 

We studied the non-equilibrium condensation process in the holographic superconductor 
by solving the time evolutions of the Einstein-Maxwell-charged scalar system in the asymp- 
totically AdS spacetime. We considered small perturbations on the Reissner-Nordstrom- 
AdS black holes, which are static solutions of this system, as the initial states. We found 
that, when the temperature is lower than the critical temperature Tc-, initial perturbations 
grow exponentially and the spacetimes settle into the hairy black holes obtained in |0. It 
is concluded that the hairy black holes are the final states for the instability of the low- 
temperature Reissner-Nordstrom-AdS black holes. We also found that the hairy black holes 
are stable against the plane-symmetric perturbations. As for the superconducting order pa- 
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rameter in the boundary theory, we clarified how it evolve in non-equilibrium process during 
the phase transition. As a byproduct, we obtained the relaxation time scale of the order 
parameter, which is inverse of the imaginary part of the fundamental quasi-normal mode. 
Finally, we studied the time evolution of the event and apparent horizons and discussed the 
their relevance to the entropy of the boundary theory. 

There are many extensions and open issues. One of them is perturbation of dynamical 
spacetimes obtained in this paper. The dynamics of transport coefficients, such as electrical 
conductivity or shear viscosity, can be known from the perturbations on the dynamical 
spacetimes. From the Dirac field perturbation, we may be able to study the Fermi surface 
properties of superconductors Il3^-p8|. These studies might deepen our understanding on 
dynamical phenomena in the AdS/CFT correspondence as well as on unknown properties 
of condensed matter system in non-equilibrium state. 

In our calculations, we imposed the plane symmetry, that is, we assumed that the space 
on (x,y)-plane is homogeneous and isotropic, as in Eq. ( p.3p . It is interesting to relax this 
symmetry. For example, we can consider the spacetime with homogeneous and anisotropic 
(x, |/)-plane as 

ds^ = — - {F{t, z)dt^ + 2dtdz) + <l'i(t, zfdx"^ + ^alt, zfdy"^ . (6.1) 

The equations of motion for this metric ansatz are given by (1 + l)-dimensional PDEs, and 
thus we can straightforwardly apply the formalism of this paper. In this setting, we can 
take into account homogeneous electric or magnetic fields along the x or ^/-directions and 
study the dynamics of the current in the boundary theory. More challenging problem is to 
take into account inhomogeneity along one direction, e.g., the ^/-direction. The equations 
of motion are given by (1 + 2)-dimensional PDEs. In the setting, we may holographically 
realize a system in which the temperature is not homogeneous and the condensed phase 
and the non-condensed phase coexist. In such an inhomogeneous system, we may observe 
the thermoelectric phenomena in superconductors |]5D|, such as Pertier effect or the Zeebeck 
effect, or the dynamics of the surface between the two phases ^,^. Another subject which 



can be studied by (1 + 2)-dimensional PDEs is the vortex solution obtained in |^2|,^: we 
must consider the ?7(l)-symmetric spacetime where the U{1) is the rotational symmetry 
around the vortex. The spacetime of the vortex solution has a ?7(l)-symmetry, which is the 
rotational symmetry around the vortex, and then the equations of motion are again given 
by (1 + 2)-dimensional PDEs. Investigation on the non-equilibrium phenomena in these 
inhomogeneous system and their comparison with theories and experiments in the field of 
the condensed matter physics may be interesting. 

Ultimately, we should study dynamics of spacetimes without any symmetry. Such a 
task may be tackled only by numerical relativity technique, which is developing rapidly 



these days pi|-p5|]. To apply this technique to asymptotically AdS spacetime, we have 



to modify the formalism for asymptotically fiat spacetime to accommodate the negative 
cosmological constant. If such a technique is established, we will be able to investigate 
on many interesting subjects concerning the non-equilibrium condensation process in the 
AdS/CFT. For example, it is expected that vortex lattice |^ are formed in a dynamical 



system from the simulation of the (non-holographic) superconducting system [|7|-[TT[]. By 
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simulating its holographic dual in the gravity side using the technique, we will be able to 
confirm further the AdS/CFT correspondence in a dynamical setup and also may clarify 
unknown properties of dynamical strongly-coupled systems. We believe that such studies 
will be fruitful for both of the quantum field theory and general relativity, and our study 
in this paper will serve as a first step toward such ambitious future issues. 
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